Optical absorption edge in one-dimensional conductors 
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The frequency-dependent conductivity is studied for both the one-dimensional Hubbard model 
and a model of spinless fermions, using a selection rule, the Bethe ansatz energy eigenstates, and 
conformal invariance. For densities where the system is metallic the absorption spectrum has two 
contributions, a Drude peak at u = separated by a pseudo-gap from a broad absorption band 
whose lower edge is characterized by a non-classical critical exponent. Our findings are expected to 
shed new light on the "far infrared puzzle" of metallic organic chain compounds. 
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Quasi-one-dimensional materials with partially filled 
electronic bands, i.e. "metals" from a naive point of 
view, usually undergo a transition to a low-temperature 
phase with a gap in the charge excitation spectrum. 
These transitions are driven by electron-electron and/or 
electron-phonon interactions. Depending on the cou- 
pling constants and the band fillings the broken symme- 
try ground states are of different kinds: spin- or charge- 
density waves, bond alternation, superconductivity. Sim- 
ple mean-field theory is able to describe these transi- 
tions, at least qualitatively. The broken symmetry in- 
duces a gap at the Fermi energy which may be displayed 
in the optical absorption spectrum. In the case of charge- 
density waves accompanied by a lattice distortion, opti- 
cal gaps have indeed been clearly observed, and beautiful 
examples have been discovered among organic molecular 
chain compounds 0], conjugated polymers j^, and inor- 
ganic chain-like materials In the case of a moving 
spin-density wave, the sliding collective mode removes 
all the finite-frequency transitions and produces a single 
peak at w = Pinning by impurities not only shifts 
this peak to a finite frequency but also restores optical 
transitions above the gap ||]. 

The mean-field scenario makes sense at very low tem- 
peratures, where interchain coupling stabilizes three- 
dimensional long-range order. In contrast, above the 
transition temperature the system may enter an essen- 
tially one-dimensional regime, which requires a more re- 
fined treatment of interactions. The real part of the op- 
tical conductivity can be written as cri(w) = 2'kD5{lj) + 
a\'^^{uj), where the frequency-dependent conductivity 
aY'^ (uj) is defined in Eq. (^) below. Giamarchi and 
Millis have used the bosonization technique and found 
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band fillings Q. Unfortunately, it is at present not clear 
to what extent their treatment is a conserving approxi- 
mation and, therefore, it is important to compare their 
results with those obtained using other techniques. Sev- 
eral authors have used exact diagonalization for the one- 



dimensional Hubbard model. In addition to finding no 
low-frequency absorption at half-filling, where the sys- 
tem is an insulator, these authors have found practically 
no low-frequency (w ^ 0) absorption away from half- 
filling, where a Drude peak at w = confirms that the 
system is metallic 0. This is not only in contrast to the 
bosonization results, but also to the two-dimensional case 
where doping away from half filling induces mid-infrared 
absorption |Q,D. Analytical calculations in the large U 
limit of the one-dimensional Hubbard model have so far 
been restricted to the electronic density n = 1 P,p^. In 
this limit the optical gap occurs at Eopt — U—At, followed 
by an absorption band extending up to ?7 -f At P^p^ . 

In this paper the optical absorption of the one- 
dimensional Hubbard model is discussed in the frame- 
work of the Bethe ansatz, which has been used previously 
by Schulz [|l^ to calculate D. Together with the optical 
sum rule, this allowed him to derive the total intensity 
of finite-frequency transitions. This intensity is generally 
small, except close to half filling for large enough (but 
not too large) values of U. We have found that these 
finite-frequency transitions are likely to be limited to a 
well-defined band above an effective optical gap, which is 
actually smallest at half filling and have also been able to 
derive the exponent for the frequency-dependence above 
the absorption edge. 

The ID Hubbard model can be written a,s H = T+UD, 
where T ~ ctHctCj+Ict + h.c] is the "kinetic en- 

ergy", D = X^j^i.T^j.i measures the number of doubly 
occupied sites, and Cj„ are electron operators of spin 

projection a at site j = 1, ...jNa, fij^^ — c^^Cja, t is the 
transfer integral, and U is the on-site Coulomb interac- 
tion. We choose a density n = N/Na in the interval 
< 71 < 1 with even N and zero magnetization. The 
Fermi momentum is given by fc^^ = 7™/2. We use units 
such that — h — 1, where — e is the electronic charge. 
The regular part of the optical conductivity can be writ- 
ten as 
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Here J 



is the current 



operator, the summation runs over energy eigenstates, 
and Lo^^o = — Eq is the excitation energy above the 
ground state |0). The formula (|l|) applies if the ground 
state is non-degenerate. This is true for an even number 
of sites Na if periodic (anti-periodic) boundary condi- 
tions are used for odd (even) values of respectively. 
The ground state is then necessarily an eigenstate of the 
parity operator, Ptt, which moves electrons from sites j 
to Na + 1 — J, j = 1, Na^ and has eigenvalues ±1. Im- 
portantly, in the present model Ptt commutes with the 
Hamiltonian and anticommutes with the current oper- 
ator. This implies immediately that the states |0) and 
J|0) have opposite parities, a key selection rule for op- 
tical transitions. The final states \v) can be charac- 
terized in terms of holons, antiholons, spinous [ p^JTs] 
and a charge-transfer band | [l6[ |. The holon and spinon 
bands describe low-energy charge and spin excitations, 
whereas the charge transfer band is associated with the 
upper Hubbard band. We use the labels a = c,s,t for 
the holon/antiholon, the spinon, and the charge-transfer 
bands, respectively and the quantum number /3 for 
distinguishing between holons {(3 — — ^) and antiholons 
{j3 = +^) 1 15 1^. In the present context, the s band 
is empty of spinous, the c band can be populated by 
holons and by zero or one antiholon, and the t band can 
have occupancy zero or one. The momentum variable 
of the different bands has the form qj — j^If, where 
/" are successive integers or half-odd integers. In con- 
trast to the case of electron bands in a periodic solid, 
the number of available momenta N* can change with 
band occupancies. These numbers are N* = iV„, N. 
N/2-Nu and N* 



ANt = 1. The former start, in 
= 0, while the latter have an onset at 
Eopt = Wt — 2e°{2kF) > 0. We show first that, as a con- 
sequence of the selection rule presented above, the tran- 
sitions of type (a) are likely to have a very small weight. 
These are excitations within the holon band which can 
be characterized in terms of single and multiple electron- 
hole excitations |lj] . Single electron-hole excitations give 
no contributions to the finite-frequency absorption, while 
multiple excitations are expected to decrease very rapidly 
in intensity. Double excitations, which would give a con- 
tribution ai{uj) (jj"^ are forbidden by the parity 
selection rule, since these transitions require symmetri- 
cal changes of holon occupancies at q and —q, leaving the 
parity unchanged. In addition, there may be higher order 
low- frequency Umklapp processes Q , but since these re- 
quire multiple excitations in the holon band, they are ex- 
pected to have extremely low spectral weight. These con- 
siderations explain why numerical studies yield so weak 
features at low frequencies l^Jl7|]. 

The weakness of type (a) transitions implies that the 
onset of type (b) transitions at Eopt represents a pseudo- 
gap which may look like a true gap in an actual exper- 
iment. Eopt, shown in Fig. 1, increases with increasing 
U and decreases with increasing density. For n = 1 it 
coincides with the Mott-Hubbard gap Emh |16|, which 
can be expressed as Emh = — 2e°(7r) and is an increas- 
ing function of U. In general, the type (b) transitions 
are also weak, except for n close to 1 and large U p3| , 
corresponding to a region of parameters where the band- 
width Wt is very small. Therefore, we have limited our 
studies to this parameter region. The important pro- 
cesses are those where in addition to a t particle a holon- 
antiholon pair is created in the c band. This generates 
an absorption band with a lower edge at w = Eopt and a 



Na-N + Nt, where Nt is the num- ^i^th ^^opt = -2[e°(0) - e°{2kF)] = 8t - 2e^. [Here 

— ^ci'^) ~ ^cC^^f), changes between ep — At for 



ber of occupied momenta in the charge-transfer band. 
Therefore, the momentum band widths are Aq^ = 27r, 
Aqs =2kF- 2n Nt/Na, and Aqt = 27r - ikp + 2t: Nt/N a- 
The numbers _i (holons), ^ (antiholons), and 

Nt are good quantum numbers which obey the sum rules 
Na-N ^ ~2Y.p=±ipNlp = N^ - 2Nt, where 

A^^ — X)/3=±i ^c,i3- The ground state is characterized 
by N'',, ^'Nt ^ and iV^_i = Na - N, with a 
symmetrical holon occupancy in the c band for momenta 
2kF < \q\ < TT. The energy bands e1{q) can be extracted 
from the Bethe-ansatz solution |l^. Interestingly, we 
find that within the parameter region of appreciable os- 
cillator strength for optical transitions [p^, the charge 



n ^ and ep 
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?(0) is 



transfer band width Wt 
very small. 

The main transitions contributing to (t['^^ (lu) are of two 
types: (a) those leaving the band filhngs N'^ i , iV'' i , 
and Nt unchanged and (b) those changing these numbers 



for 71 ^ 1.] The use of the quasi- 
particle/electron representation of Ref. ||l^ reveals that 
these transitions are associated with one zero-momentum 
electron - hole process. Moreover, creating a t parti- 
cle means enhancing double occupancy by 1, at least 
in the large U limit. Since the current operator can- 
not produce more than one doubly occupied site, there 
will be almost no spectral weight above the upper edge 
at = Emh + 8t + « Emh + 8t, illustrated by the 
dashed line in Fig.l. 

The spectral range of optical absorption deduced above 
on the basis of the holon and charge transfer excita- 
tion spectra fully agrees with large U expansions |l^] 
and numerical calculations [[7|JTl|]. A much more diffi- 
cult problem is the spectral line shape. Here, we limit 
ourselves to the region immediately above threshold, 
using the hypothesis of conformal invariance. In this 
region the transition energies have a finite-size depen- 
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dence - Eq - Eopt = Ea,.=±i + Oil/Na) 

where AJ^ are U and n dependent positive parameters. 
The momentum is the same as in the ground state, i.e. 



P. 



N„ l^c 



= ±1 



lA'-^ = 0. By analogy to the case of 



the two- velocity spectrum of Ref. |18|] , we interpret this 
as a low-energy {uj — Eopt) three- velocity conformal- field 
theory, whose vacuum is |0} with energy measured from 



En+E. 
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We have verified that the exponents for the lem 



nent from the data. In the past there has been consider- 
able discussion concerning the phenomenon of suppressed 
far infrared conductivity observed in organic chain com- 
pounds, despite their large dc conductivity | |2l[ |. Very 
recently, the optical absorption spectrum of large sin- 
gle crystals of Bechgaard salts has been measured over a 
huge frequency range, confirming the far infrared prob- 
22 1 . A Drude peak corresponding to a very small 



conductivity do not depend on whether they are calcu- 
lated for vt — >• or for vt = 0- The physical fields 
are of the general form (p{x,t)e'^^°''** . The anomalous 
dimensions of the (j){x, t) fields involve the parameters, 

^aa' = '^".a' +J2l=±ll'^aa'{qFaJqFa') (with Qps = kp 

and qpt — 0). Here j = 0, 1 and the two-particle forward- 
scattering phase shifts, ^aa' iq, q'), can be extracted from 
the Bethe-ansatz solution The critical behavior of 

the conductivity (|l|) can then be calculated, with the re- 
sult a\^^ {lo) ~ c{uj — Eopt\ at the onset and a critical 
exponent 
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The exponent C: shown in Fig. 2, approaches ^ as 
[/ — > 00, and for densities n < 1 it increases with decreas- 
ing values of U . For n ^ 1 we find C = \ ^'^d Wt = for 
all finite values of U . Putting = and n = 1 in the 
X and t dependent current - current (or charge - charge) 
correlation function and taking the Fourier transform, 
we again obtain C = ^ at the threshold of optical absorp- 
tion. Therefore, we believe that at n = 1 the expression 

(7^^ {ijj) = C^w — Emh^ is valid for all finite values of 

U . For the particular case of n = 1 the same method 
allows us to evaluate the leading frequency dependence 
close to the upper edge of the absorption band. We find 



carrier density is separated by a gap of about 0.01 eV 
from an absorption band of a width of about 0.6 eV. The 
onset of absorption above the gap is very sharp at low 
temperatures. In the framework of the one-dimensional 
Hubbard model these results can be understood if the 
system has a nearly half-filled band, in which case the op- 
tical threshold is essentially given by the Mott-Hubbard 
gap. The parameters t « 0.25ey and U /t w 1.5 allow to 
reproduce the gross features of the experimental spectra. 

While the simple Hubbard model may be sufficient for 
reproducing the overall absorption spectrum of some of 
the one-dimensional conductors, additional terms such 
as the Coulomb interaction between nearest neighbor 
sites would have to be considered for a more detailed 
description of the spectrum, for instance for the exci- 
tonic enhancement at the absorption edge p^ . Since 
the extended Hubbard model is not solvable by the 
Bethe ansatz, we have considered the model Hamiltonian 
H — —t^j[c]cj+i+h.c]+V^^ fijUj+i, where c] creates 

aspinless fermion at site j — 1, Na, fij — c^Cj, and V is 
the nearest-neighbor Coulomb interaction. For a density 
n ~ — this model has a metal - insulator transition 
for V > 2t. As for the Hubbard model, in the region 
V > 2t we find almost no spectral weight in ai{u!) for 
Lu < Eopt, where Eopt is a ^-dependent optical gap. We 
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again obtain a critical behavior cr['^^(w) = c(^uj — E, 
(w) = c(Emh + 8t-u;y. This symmetry between low values of {uj - Eopt), with ^ 5 both as n 



the lower and upper absorption edges for n = 1 was also 
found in Refs. jl^,|lj]. Moreover, expanding the n ~ \ 
expression (28) of Ref. |l0| around the edge energies leads 
precisely to the same uj power-law dependences with ex- 
ponent \. We notice that the constant C has to vanish 
for all densities as C/ — s- or [/ ^00, since in these limits 
the Drude peak exhausts the conductivity sum rule [p^ . 

It is worthwhile to compare our predictions with 
experiments on one-dimensional conductors. Unfortu- 
nately, as a general rule these materials cannot be rep- 
resented in terms of the Hubbard model alone, but 
additional interactions have to be taken into account. 
An exception seems to be the organic chain compound 
{TTM - TTP)l3, which displays the properties of the 
one-dimensional Hubbard model at half filling with a rel- 
atively small value of [/ (C/ w 3.5i): an optical gap of the 
expected size |19 and a Pauli-like susceptibility |2^]. The 
optical absorption is consistent with a smooth onset at 
threshold, but it seems difficult to extract a critical expo- 
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and V > 2t and for all densities as F — + 00. For n = ^ 
and large V the corresponding absorption extends from 



UJ 



V - At to UJ ^ V + At. 



In summary, we have discussed the frequency- 
dependent conductivity of the one-dimensional Hubbard 
model (and of a model of spinless fermions with near- 
est neighbor interaction), using a parity selection rule 
and known properties of the exact eigenstates. In the 
insulator phases the optical absorption has a square-root 
behavior at the onset. For n ^ 1 (n 7^ 1/2 and V > 2t) 
the spectrum consists of a Drude peak ai u — and an 
absorption band above a pseudo-gap, which increases as 
the system is doped away from the insulating state. The 
behavior at the onset of the finite-frequency transitions 
is described by a critical exponent that depends both on 
the interaction strength and the density. The intensity of 
the transition across the gap is appreciable only when the 
charge-transfer bandwidth is very narrow. This is remi- 
niscent of the X-ray edge singularity which depends very 
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sensitively on the recoil of the hole Our results are 
expected to shed new light on the old problem of optical 
absorption in one-dimensional conductors. 

We thank D.K. Campbell, P. Horsch, and J.M.B. 
Lopes dos Santos for illuminating discussions. 
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FIG. 1. The gap Eopt, (a) as function of U for different 
values of n and (b) as function of n for different values of U. 



FIG. 2. The exponent (a) as function of U for different 
values of n and (b) as function of n for different values of U. 
The straight lines ai ( — 1/2 refer to the (a) n — > 1 and (b) 
17 ^ 00 limits. 
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